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We propose a method of measurement of extragalactic magnetic fields in observations of TeV 
7-rays from distant sources. Multi-TeV 7-rays from these sources interact with the infrared photon 
background producing secondary electrons and positrons, which can be just slightly deflected by 
extragalactic magnetic fields before they emit secondary 7 rays via inverse Compton scattering of 
cosmic microwave background photons. Secondary 7-rays emitted toward an observer on the Earth 
can be detected as extended emission around initially point source. Energy dependent angular profile 
of extended emission is determined by the characteristics of extragalactic magnetic field along the 
line of sight. Small magnetic fields B < 10"^'^ G in the voids of the large scale structure can be 
measured in this way. 

PACS numbers: PACS: 95.85.Pw, 98. 62. En 



The problem of the origin of 1 — 10 /iG magnetic fields 
in galaxies and galaxy clusters is one of the long standing 
problems of astrophysics/cosmology (see [l| for a review). 
Such fields are thought to be produced either via dynamo 
mechanism [3| or via compression of primordial magnetic 
field during the large scale structure (LSS) formation [l|. 
In either of these two scenarios one assumes the existence 
of a small " seed" primordial magnetic field of cosmologi- 
cal origin. The assumed strength of the primordial mag- 
netic field and its correlation length are strongly model 
dependent. Moderate experimental limit on the present 
day strength of primordial magnetic field B < 10~®G 
comes from the limit on rotation measure of emission 
from distant quasars, assurning (largely uncertain) cor- 
relation length Ic ^ 1 Mpc [3]. Similar restriction comes 
from the analysis of anisotropics of the cosmic microwave 
background (CMB) ^]. 

Recently a significant progress was achieved in calcu- 
lations of the tree-dimensional structure of extragalactic 
magnetic fields (EGMF) [5|, |6| . These calculations are 
based on the numerical simulations of LSS formation. 
The main result of the magnetic field simulations is that, 
similarly to the LSS, the magnetic field structure is fil- 
amentary, with large voids, filaments and knots. The 
strength of the primordial magnetic field in these simu- 
lations is chosen in such a way as to reproduce the mea- 
sured fields B ^ 1 nG in Galaxy clusters. New simula- 
tions allow to estimate the cumulative volume filling fac- 
tor V{B) for EGMF with the strength close to a typical 
value B. One should note that the results of simulations 
by different groups differ dramatically. In particular, in 
Ref. [1 it is claimed that V(S > 10^^ G) ~ 0.15, while 
V{B < 10-12 G) ~ 0.07. To the contrary, in Ref. [6] 



V{B > 10-9 G) ~ 2 X 10-4 while V{B < lO-^^ G) ~ 0.7. 
Recent independent calculation of ref. ^TJ suggests that 
B > 10-9 G fields fill ^ 1% of volume, which is somewhat 
in between of the results of Ref. [5| and Ref. [g| . 

Obviously, the best way to resolve the above contro- 
versy of numerical calculations is to measure the EGMF 
outside the galaxies and galaxy clusters. In this paper we 
propose a way of measurement of EGMF which is sensi- 
tive to the very weak magnetic fields B < 10- ^^ G and, 
therefore, could be used to probe, for the first time, the 
magnetic fields in the voids of the LSS. 

The idea is to look for the extended TeV 7-ray emis- 
sion around distant sources. It is known that multi-TeV 
7-rays from these sources could not reach the Earth di- 
rectly because of the pair production on the extragalactic 
(infrared) background light (EBL). The pair production 
on EBL should lead to the exponential suppression of the 
7-ray flux from the source, 



F{E,„) ^ Fo{E^„)exp[~T{E,„,z)] 



(1) 



Here F{E^g) is the detected spectrum, Fo{E-yg) is the 
initial spectrum of the source and T{Ejg,z) is the optical 
depth with respect to the pair production on EBL, which 
is a function of the primary photon energy E^^ and of the 
redshift of the source z [8| 



r = 5^10 



z 

or 



PIR/O 



(5 ^ 10) nW/(ni sr) 



E. 



70 



10 TeV 



(2) 



(we allow for an uncertainty of about a factor of 2 in the 
density of the infrared/optical background, Pir/o, taking 
into account the recent controversy in the measurements 
of COBE/DIRBE and the upper limits on pinjo imposed 
by the TeV observations of distant blazars [9| ) . 
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FIG. 1: Schematic representation of formation of extended 
TeV emission around a distant point source due to the cascade 
on the EBL in small EGMF. 



pairs of the energy Ee produced in inter- 
actions of multi-TeV 7-rays with EBL photons produce 
secondary 7-rays via inverse Compton scattering (ICS) 
of the CMB photons to the energies 



4 Ei 

Ej = -ecMB — % — 1-2 



E. 



70 



40 TeV 



TeV 



(3) 



(Here ccmb = 6 x 10"'' eV is the typical energy of 
CMB photons. We also have assumed that the energy 
of primary 7-ray is E-y^ ~ 2Ee)- Upscattering of the 
infared/optical backgound photons gives sub-dominant 
contribution to the ICS because the energy density of 
CMB, PC MB '^ 10 -^ 100p/ii;/o, is much higher than the 
density of the infrared/optical background. 

It is conventionally assumed that as soon as trajec- 
tories of electrons are deflected by the EGMF, the sec- 
ondary photons are not emitted in the direction toward 
the Earth and do not reach the telescope. However, the 
above way of reasoning is correct only if one calculates 
only the reduction of the flux of photons emitted initially 
exactly in the direction of observer. Deflections of elec- 
trons produced by the 7-rays which were initially emitted 
slightly away from the observer, could lead to "redirec- 
tion" of the secondary cascade photons toward the ob- 
server. This effect leads to the appearance of additional 
"cascade" contribution to the source flux. 

The cascade contribution can be distinguished from 
the " primary" source contribution with the help of imag- 
ing capabilities of the TeV Cherenkov telescopes. The 
secondary cascade 7-rays are emitted at large distances 
from the source, from an extended conical region (blue 
shaded region in Fig. [1]). The opening angle of the cone 
is determined by the deflection angle of secondary elec- 
trons, (5, its vertex coincides with the point source, and 
its height is equal to the mean free path D^^ — D /t ~ 
10 -=- 20[^T,o/40 TeV]"^ Mpc {D is the distance to the 
source) of the multi-TeV photons through the EBL. Since 
the line of sight is located inside of the cone, the cone 
appears as an extended emission around the initial point 
source (Fig. [T]). 

The mean free path of the secondary electrons with re- 
spect to the ICS of CMB photons, Ae = {(TTncMB)~^ — 
1 kpc ((TT = 6.7 X lO^^^cm^ is the Thompson cross- 
section and ncMB — 411/cm^ is CMB photon density), 
is much shorter than D^„ and one can assume that elec- 



trons immediately re-emit 7-rays via the ICS. The de- 
flection angle of electrons, i5, is determined by the ratio 
of Ag to the Larmor radius in the EGMF of the strength 
B, tl ^ 6.6 X 10^2 [Se/20 TeV] [lO'^^ G/B] cm. As- 
suming Ae ^ tl one can find 5 ~ ^b/tl — 
3.2° [20 TeV/Se] [S/IO^^^ qj a simple geometrical 
calculation shows that the angular size of extended emis- 
sion produced via redirection of cascade photons toward 
the Earth is (see Fig. [1]) 
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(4) 



(throughout the paper we use the relation between the 
primary photon energy, E^^, the energy of secondary 
electron, E^ and of the secondary photon, E^ given by 
Eq. USD). 

Comparing the distance covered by the photons which 
come directly from the point source with the distance 
covered by (primary and secondary) cascade photons, 
one can find (from the same geometrical calculation) 
the typical time delay of the secondary photons, td — 
S^D{t — 1)/2t^. Contrary to the point source fiux, the 
flux in the cascade component is not variable, since it is 
equal to the point source flux averaged over a large time 
scale td (e.g. for ,5 ~ 5°, z = 0.03 and E^„ - 40 TeV, one 
flnds td ^ 10^ yr). 

Since each secondary cascade photon carries only a 
small fraction of the primary photon energy (see Eq. 
([3])), the flux of the extended source, integrated over 
the region of the angular size 9ext is a small fraction 
a ~ 2E~jE-y^^ ~ 0.06 [E-y^/AQ TeV] of initial point source 
flux {Fo{Ejg)) (angle brackets signify averaging of the 
point source flux over the time scale td). The extended 
source flux is suppressed at high energies because the 
secondary cascade photons can themselves produce pairs 
in interactions with EBL photons. The suppression fac- 
tor is exp{—T{Ej,z)) where r is given by ([2]). Apart 
from the absorption at high energies, the spectrum of 
the extended source is modified if only a fraction of the 
mean free path of the primary 7-ray, Dj^ , goes through 
the voids with small magnetic field B. This fraction 
can be expressed through the volume filling factor of 
V{B), as Dvoids/D-yg ~ V^^^. Combining all the cor- 
rection/suppression factors one can find a relation be- 
tween the average absorbed point source flux {F{Ejg)) 
and the flux of extended source integrated over the re- 
gion 9 < 9ext- 
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Cascade electrons loose about half of their en- 
ergy after upscattering of many CMB photons over 
the energy attenuation distance Ae ^ (Ee/E^)\e ~ 
20 [20 TeV/E'e] Ae- ICS emission produced by these elec- 
trons over the distance Ae > Ae forms a powerlaw tail 
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FIG. 2: Surface brightness profile of Extended source produce 
by the cascade on EBL. 



around the core of the extended source, 9 < 9ext, as it is 
shown in Fig. [2l Since for electron propagtion distances 
d < Ae the deflection angle scales linearly with d, the 
surface brightness profile of the tail can be found by in- 
tegrating the ICS energy loss along electron trajectory, 
S{d) ~ ^"^(l + 29/&ext), where Qext is the opening an- 
gle of a cone into which the ICS photons produced by the 
cascade electrons over the energy attenuation length are 
emitted, A = AJvl ^ 6° [20 TeY/E^f [B/IQ-^^ G] . 
Using Fig. [T]one finds. 
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(6) 



The fraction of the total flux of extended source inte- 
grated over the region 9 < Oext is about the initial point 
source flux {Fo{E^g)) modified at high energies by the 
effect of absorption of the secondary 7-rays and by the 
volume-filling-dependent factor V^^^, 



FextiEy, Qext) — V ' 



f^riE^g 



1 



r{E^,z) 



{F{E^,)) (7) 



(compare with Eq. ^). 

Up to now we have assumed that the primary 7-rays 
are emitted by the point source isotropically. However, 
most of the distant sources of TeV 7-rays are blazars. 
The primary 7-ray emission from blazars is known to be 
concentrated along the direction of the jet with a typ- 



ical opening angle Q 
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10 is the 



typical bulk Lorentz factor of the jet). There is no differ- 
ence between isotropically and anisotropically emitting 
source till the development of the cascade on EBL leads 
to redistribution of directions of secondary 7-rays within 
a cone with an opening angle smaller than the opent- 
ing angle of the jet (see Fig. [1]). However, as soon as 
the deflected cascade photons are emitted into a cone 
with opening angle larger than Qjet , the flux of extended 
source becomes suppressed. This leads to a steepening 
of the surface brighness profile of the extended source at 



the angles 9 > 9cut where 



D~D^ 



-e 



e 



jet 



jet 



1 



(8) 



Suppose that extended sources with characteristic 
(energy-dependent) surface brighness profiles shown in 
Fig. [2] are found around a number of distant point 
sources. In this case measuring the extended source pa- 
rameters one can constrain the characteristics of EGMF 
along the source directions. For example, combining Eqs. 
(HI) and (O one finds that measurements of 9ext{E^) and 
FextiEj, 9ext) provides a measurement of the volume fill- 
ing factor V{B) of EGMFs of particular strength B: 



B~ 10-'-'t{E^„,z) 
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(9) 



Othervise, from ([6]), ([7|) one finds that information about 
B, V{B) can be extracted from the measurement of 

QextyE^), FextyE^, Qext)'- 
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(10) 



Determination of characteristics of EGMF from either ^ 
or PH)) constitutes the essence of the proposed method 
of measurement of EGMF. 

Maximal EGMF strength which can be probed is found 
from observation that 9ext should not be larger than 9cut 
(to be measurable). Since 9ext ~ E~^/'^ (see (|3])), the 
maximal EGMF strength is determined by the condition 
that 9f,xt{Emax) < 9 cut where E^ax is the highest energy 
at which a source is detected. Combining ([4]) and ([8]) one 
finds 



Br, 



3x10" 



r(B^o , z) 



t{Ej„,z)-1 



E„ 



ii/2 r 



10 TeV 



e 



jet 



G, 

(11) 
The weakest magnetic fields which can be probed is 

found from observation that Qext can not be measured 
if it becomes smaller than the size of the point spread 
function of the telescope, 9psF- Since Qext ^ E~^ (see 
([6])), the largest source extension is achieved at lowest 
energies. Taking into account that for present generation 
instruments the low energy threshold is situated at the 
energies E,nin ~ 100 GeV and that the typical point 
spread function of a Cherenkov telescope has the size 
9psF --0.1°, one finds from dH) and ^ 
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As an example of implementation of the proposed 
method of measurement of EGMF, let us find the con- 
straints on V{B) imposed by the non-observation of ex- 
tended emission around Mkn 501 by HEGRA. An upper 
limit on the extended emission flux within 0.5° around 
Mkn 501, Fext,o.5° < 10~^ Ffiare {Ffiare IS the flux from 
Mkn 501 in the flaring state) was derived by HEGRA 
collaboration [lO|). Unfortunately, neither the energy de- 
pendence of the derived upper limit nor the assumed sur- 
face brightness proflle are reported in [lO| . We assume in 
the following that the restriction was put on the flux in 
the core of the extended source, 6 < 9ext- The quiescent 
flux of Mkn 501 is about 10 times lower than the flux 
in the flaring state which means that HEGRA observa- 
tion imposes a restriction on the extended emission flux 
at the level of Fe:.t(0.5 TeY^Oext) < 0.01 (F(0.5 TcV)) 
The primary 7 rays which produce extended emission 
at Ej ~ 0.5 TeV have energies Ejg ~ 26 TeV. Optical 
depth for such 7 rays is large, r(26 TeV, 0.03) =4 — 8 
(see ([2])) even for the relatively nearby source Mkn 501 
{z — 0.03). The spectrum of Mkn 501 is characterized by 
the photon index F = 2.3 ± 0.2 [ll|, which means that 
(Fo(26 TeV)) ~ 0.3 (i^o(0.5 TeV)) or 



Fext{E^,eext) < Om {F{E^,,)) e 



•(Sto'^) 



(13) 



Using Eq. ([9]) one flnds that the magnetic field 
strength probed by the study of extended emission at 



0.5° is 
dH) one 



Ej ~ 0.5 TeV at the angular scale 6ext — 
S - (3 ^ 6) X 10-13 Q Substituting ^ into 
finds 

V [B ~ (3 ^ 6) X 10-1=^ G] < 0.4 (14) 

which is much above the predictions of Ref. [5J {V{B ~ 
10""'^^ G) ^ 0.07) and is just about the prediction of Ref. 
{V{B ~ 10-12 Q) ^ 0.7). The above example shows 
that with modern telescopes, like HESS, MAGIC and 
VERITAS, one has a real possibility to detect or put 
tight constraints on the EGMF in the 10-^^ - lO-^^ G 
range if one makes a systematic search of extended emis- 
sion around a large number of point sources at different 
energies and different angular scales. 

Electromagnetic cascades of multi-TeV photons on 
EBL were considered before in a different context. If 
the cascade happens directly in or near the source, sec- 
ondary electrons and positrons will be completely ran- 
domized in relatively large magnetic fields, B > IQ-^ G, 
producing R > 1 Mpc-size halo around the source [12| . 
In an opposite case, when EGMF is extremely small, 
B < 10-1* G, all the cascade will proceed in the for- 
ward direction and just contribute to the point source 
flux (even during the flares) [13|, [ij] . Contrary to both 



those cases we are interested here in magnetic fields in 
the intermediate range 10-^^ G< B < IQ-^^ G, in the 
voids of LSS, and stress that such fields can be measured 
by detection of the time independent extended emission 
structures around observed TeV point sources. 

To summarize, we have proposed a method of mea- 
surement of extra-galactic magnetic fields with the help 
of TeV Cherenkov telescopes. The idea is to look 
for extended emission produced by cascading of multi- 
TeV 7-rays emitted by distant point sources. The ex- 
tended emission produced by the cascade in relatively 
weak EGMF is expected to have a characteristic energy- 
dependent surface brighness profile. Measuring the pa- 
rameters of this profile enables one to determine the 
characteristic of EGMF from Eqs. jS]) and/or (fTO| . 
The method is sensitive to EGMF strength in the range 
10-16 G<B < 10-12 G. 
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to this manuscript. 
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